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Abstract: Though the operator product expansion is applicable in the calculation of 
current correlation functions in the Euclidean region, when approaching the Minkowskian 
domain, violations of quark-hadron duality are expected to occur, due to the presence of 
bound-state or resonance poles. In QCD finite-energy sum rules, contour integrals in the 
complex energy plane down to the Minkowskian axis have to be performed, and thus the 
question arises what the impact of duality violations may be. The structure and possible 
relevance of duality violations is investigated on the basis of two models: the Coulomb 
system and a model for light-quark correlators which has already been studied previously. 
As might yet be naively expected, duality violations are in some sense "maximal" for zero- 
width bound states and they become weaker for broader resonances whose poles lie further 
away from the physical axis. Furthermore, to a certain extent, they can be suppressed by 
choosing appropriate weight functions in the finite-energy sum rules. A simplified Ansatz 
for including effects of duality violations in phenomenological QCD sum rule analyses is 
discussed as well. 
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1 Introduction 

In phenomenological applications of QCD, the operator product expansion (OPE) [1, 2] 
plays an essential role. It is for example employed in applications of QCD sum rules 
[3] or the analysis of the r hadronic width [4]. Considering in particular two-point cor- 
relation functions of mesonic currents, strictly speaking the OPE is only appropriate in 
the Euclidean domain where s = q 2 is sufficiently large and negative, with q being the 
four-momentum transfer between the two currents. In the physical, Minkowskian domain 
with positive s, generally bound states or resonances are present which are governed by 
long-distance physics and which cannot be captured by a short-distance expansion. 

In the Euclidean domain, a duality exists between descriptions of the system in terms 
of the quark and gluon degrees of freedom of QCD, or the physical hadronic degrees of 
freedom. Upon analytic continuation of correlation functions towards the Minkowskian 
region, due to the presence of bound-state poles, violations of the quark-hadron duality 
are expected to show up [5-8]. 1 Since until now no analytic solution to QCD has been 
found - and it is unlikely that this will happen anytime soon - the precise functional form 
of duality violations (DVs) remains unknown. This necessitates to resort to the study of 
models in order to investigate their influence in phenomenological analyses of QCD. This 
is particularly pressing when employing so-called QCD finite-energy sum rules [10, 11], 
because here contour integrals in the complex energy plane have to be computed down to 
the Minkowskian axis. 

Based on the models presented in refs. [5, 6, 8, 9], the influence of duality violations 
on hadronic decays of the r lepton has been investigated in recent years [12-16], because 
they may play a role in precision determinations of fundamental QCD parameters, like 
the strong coupling a s [17-21]. For example it has been observed that the compatibility 
of values for the gluon condensate {a s GG), extracted from fits to vector and axial vector 
r decay spectra [17, 20], is not very satisfactory. This poses the legitimate question if 

1 For a review the reader is also referred to ref. [9]. 
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the inclusion of DVs in a consistent analysis may improve the situation regarding the 
compatibility of condensate parameters extracted from different channels, and in how far 
this would affect the resulting value of a s . 

Though the particular model for vector and axialvector correlation functions employed 
in the analyses [12, 13] is based on features of large- N c QCD [22] and Regge theory [23, 24], 
it is not directly rooted in fundamental QCD. Thus, it appears important to detect other 
systems which allow for an independent investigation of DVs, in particular their functional 
form as well as their possible impact on phenomenological analyses. A system that allows 
for an analytical treatment, and bears close relation to quarkonia in QCD in the limit of a 
heavy quark mass, is the Coulomb system. In the context of QCD moment sum rules, to 
some extent this example was already studied in ref. [25, 26]. Furthermore, in relation with 
hadronic decays of B mesons, duality-violating contributions resulting from charmonium 
states were also considered in refs. [27, 28]. 

In the first half of this article, it will be explored in more detail what can be learned 
from the Coulomb system regarding violations of duality. To this end, the notion of duality 
will be adopted in a more general context. The point of view taken in the following will 
be that a certain expansion is employed which in general can only be considered to be of 
an asymptotic nature. In our case this refers to the OPE, and duality is supposed to imply 
that the asymptotic expansion provides an acceptable representation of the full function. 
On the other hand, in some regions of the expansion, formally exponentially suppressed 
contributions may become relevant, and those will be identified with DVs. 2 

One finding of the study of the Coulomb system below will be that due to the presence 
of zero-width bound-state poles on the real s-axis, the effects of DVs turn out to be par- 
ticularly strong. The same conclusion has been drawn on the basis of the 't Hooft model 
[29]. For this reason, the example of the Coulomb system will mainly be used to gain 
further general insights into the structure of DVs. For a more quantitative analysis, in the 
second half of this work recourse will again be taken to a simplified version of the model 
employed in refs. [12, 13] in the case of a finite width for the resonances. The finite width 
entails that the bound state poles move to the unphysical region of the complex s-plane 
and consequently, DVs in the physical region turn out to be less pronounced. 

In refs. [13, 14] a simplified Ansatz was advocated in order to incorporate DV contribu- 
tions into phenomenological QCD analyses. This Ansatz is supposed to be admissible if the 
resonances are sufficiently broad and thus lie sufficiently far away from the Minkowskian, 
physical axis. Towards the end of this article, arguments will be given why a related Ansatz 
should also be considered on an equal footing which is applicable in part of the complex 
energy plane, and the corresponding Ansatz is provided. 

2 Sum rules for the Coulomb system 

In the following, sum rules for the quantum mechanical Coulomb system will be set up, 
which are analogous to the sum rules studied in QCD [3]. Quantum mechanical sum rules 

2 This somewhat broader perspective could also be applied to the perturbative expansion with the "du- 
ality violations" being the power corrections due to the presence of QCD vacuum condensates. 
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had already been investigated in the early years of QCD sum rules [30-33], as analytical 
examples in which the techniques applied in QCD sum rule analysis could be studied and 
tested. Here, this route shall be followed for the question of duality violations. 

Before going into the particular example of the Coulomb sum rule, let us set up the 
framework of quantum mechanical sum rules in more general terms [33] . Consider a particle 
of mass m under the influence of a potential V{x). The Schrodinger equation for stationary 
states takes the form 



Hif>(x) 



ip(x) = Eif)(x). (2.1) 



Let ip a {x) be the eigenfunction corresponding to the eigenvalue E a . The resolvent operator 
G(z) of H is defined as 

G{z) = \h-zi] , (2.2) 

where z is an arbitrary complex number. Its matrix elements in the position representation 
are given by 

C (i , y; „ = m , M - £ MM®. + / d£ , ( ,3) 

where the sum runs over the discrete spectrum, the integral is taken over the continu- 
ous spectrum and Pcont 

(x, y; E) is the spectral density corresponding to the continuous 

spectrum. 

Let us introduce the function D{E) which plays the analogous role of a physical cor- 
relation function in QCD sum rules, for example the Adler function [34]: 



±G(x,y;E) 



(2.4) 



x=y=0 



D(E) : 

Using eq. (2.3), D(E) can be expressed as 

n(F\ = X^ l^(°)l 2 f gcogtCgQdg f P (E')dE' 

[ ' ^ (E a — E — *0) 2 J (E'-E- iO) 2 J (E'-E- iO) 2 ' 1 ' 

with appropriate limits of integration, and where we furthermore have introduced the full 
spectral function 

Pcont 

(E) = \M0)\ 2 S(E - E a ) + Pcont (E) , (2.6) 

a 

containing both the discrete as well as the continuous spectrum. 

Let us now particularise the general expressions to the Coulomb problem. The ensuing 
sum rules turn out to be closely related to heavy-quark sum rules in QCD. Consider two 
particles of mass m coupled by a Coulomb potential of the form 

V(r) = - " . (2.7) 



-3- 



The Coulomb potential is written with a general coupling a, though one may well regard 
this as the leading term in the heavy-quark potential containing the QCD coupling a s . 
Solving the relevant Schrddinger equation, the radial Green function for the Coulomb 
potential is found to be [25] 3 

G(r, 0; E) = ^ e~ kr T(l - A) U{1 - A, 2; 2kr) , (2.8) 

with 

am 



A = — and k = yj- m{E + iO) . (2.9) 

U(a, (5\ z) is the confluent hypergeometric function. G(r, 0; E) is singular in the limit r — > 0, 
but the physical correlation function D(E) remains finite and is found to be 



d 



dE 



mX 



r=o 4-rra 



D(E) = -f-G(r,0;E) = 1 + 2A + 2AV(1 - A) , (2.10) 



with ip(z) = d[lnr(z)]/dz being the logarithmic derivative of the Gamma function. The 
function D(E) has poles at A n = n, n = 1, 2, . . . which translate to the bound-state energies 
E n of the Coulomb system with 

„ a 2 m 

Expanding D(E) around the position of the bound states, E = E n — e, from the leading 
1/e 2 singularity the square of the wave function at the origin can be extracted, 

l<M0)| 2 = (2.12) 

which of course agrees with well known results in textbook quantum mechanics [36]. 

Besides the discrete spectrum for E < 0, in the case of the Coulomb problem also 
a continuous spectrum is present for E > 0. Taking the imaginary part of eq. (2.3), 
one observes that the spectral density p(E) is related to the imaginary part of G(r,0;E), 
namely 

p(E) = -JmG(0,0;E + i0). (2.13) 

Again, also Im G(r, 0; E) is finite in the limit r — > since it is a physical quantity. Calcu- 
lating p cont {E) from eq. (2.10), yields 

2 

oim 

Pcont(E) = — = ==-, E > 0, (2.14) 



4tt 



which is known as the so-called Sommerfeld factor [37]. Putting everything together, 
a finite-energy sum rule, analogous to the case of hadronic r decays, for the quantum 



3 As a matter of principle, also higher-order QCD corrections could be included in a systematic way. See 
e.g. ref. [35]. However, for the present purposes, it suffices to stay at the leading order. 
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mechanical Coulomb problem can be written down: 

E 



R W (E ) = J w{E)p(E)dE = j> w(E) G(0, 0; E) dE 



Ei E 



= j>W{E)D(E)dE , (2.15) 

Eo 

which defines the moment R w (Eq) corresponding to an arbitrary analytic weight function 
w(E), and the weight function W(E) is given by 

E 

W(E) = j w(E') dE' . (2.16) 

E 

The central idea of finite-energy sum rules is to employ a phenomenological representation 
of the spectral function p(E) on the left-hand side of eq. (2.15) and a theoretically moti- 
vated, like the OPE, for D(E) in the contour integral on the right-hand side. Equating 
both sides then allows to infer information on the physical spectrum from the theoretical 
description, or to extract theoretical parameters from the phenomenological information. 
As most probably in practical applications the theoretical expansion of D(E) on the right- 
hand side is only of an asymptotic nature, the consequences for the finite-energy sum rule 
need to be investigated. 

3 Asymptotic expansion and numerical analysis 

For sufficiently large positive or negative energy E, that is |A| < 1, the ^//-function appear- 
ing in D(E) of eq. (2.10) has a convergent expansion [38]. Therefore, in this energy region 
no duality violations are expected to occur. They may, however, arise in energy regions 
where tp'(z) only has an asymptotic expansion, which is the case for \z\ = |1 — A| going to 
infinity, while | argz| < tt. This happens if the energy is close to the continuum threshold 
E = 0. In this region, the asymptotic expansion of ip'(z) is given by 



/ / \ 1 1 \ — ^ / i i \ / \ 

~ z + 2^2 +Z^^TT (*->ooin|argz| < vr) , (3.1) 

n=l 

where -E?2n are the Bernoulli numbers. Even though strictly speaking the asymptotic ex- 
pansion (3.1) should be valid in the full region | argz| < tt, due to the poles of ip(z) for 
z = 0, — 1, —2, . . ., the asymptotic expansion becomes very inefficient if z approaches these 
poles. Quite generally, in the left-half imaginary z-plane, that is Hez < 0, the asymptotic 
expansion can be greatly improved by applying the so-called reflection relation [38] 



sin (ttz) 



z-^oo 1,1, \ - ti^n K 



T + ^ + (*»«>). (3-2) 

n=l x ' 



- 5 - 



A more quantitative account on the improvement achieved by including the additional term 
in the expansion of (3.2) will be given in the second example of the next section. 

As the asymptotic expansions of ip'( z ) an d —ip'(—z) + 1/z 2 are identical, it is clear 
that the additional term tt 2 / sin 2 (7rz) has to be exponentially suppressed. This can easily 
be verified by rewriting the sin-function in terms of exponentials, which yields 



sin (ttz) 



47r 2 e ±2niz + u±**iz\ for Im z > o (3.3) 



On the other hand, close to the poles of ip'(z) the exponentially suppressed term provides an 
essential contribution as it gets enhanced by the poles contained in 1/ sin 2 (7rz) for Im z = 0. 
Thus, in the general terminology employed in the present work, the additional term (3.3) 
can be considered a duality-violating contribution. 4 

In the following numerical analysis, the FESR (2.15) shall be investigated in a way 
such that the asymptotic expansion can be performed in a most transparent fashion. A 
convenient choice of the complex integration contour to this end are paths with constant 
|1 — A| = |z| = zo- The corresponding energy E z {tp) is then given by 

^ , \ a 2 m , 

^> = - 4(V-D" (3 ' 4) 

with the parametrisation z = ZQe lLp . Because of the square in the denominator, in the 
full range of < ip < 2ir, the contour E z ((p) covers the complex Emplane twice. However, 
only the range po < p < 2tt — ipo, where the angle <po = arccos(l/2;o), and which contains 
the section on which the DV term (3.3) should be included in the asymptotic expansion 
of ijj'(z), provides the correct solution to the corresponding transcendental equation. Of 
particular interest are the two points 

r-, ^ a 2 m —Ei 

E ^ E ^)= 4(?=i) = (^i)' (3 ' 5) 

a 2 m Ei 

where E z ((p) is real, and which have also been expressed in terms of the lowest-lying bound- 
state energy E\. Eq. (3.6) shows that in order for the contour not to hit a bound-state 
pole, zq should not be integer, and that depending on zq, poles with n > int(zo) + 2, where 
int(zo) denotes the integer part of zq, need to be included on the phenomenological side of 
the sum rule, as they lie inside of the integration region. 

As a particular example, the sum rule (2.15) will now be investigated numerically for 
the trivial weight function w(E) = 1 and the set of parameters a = 1/2 and m = 10. 
This leads to the lowest ground state being at Ei = —5/8. The integration contour in 



4 The appearance of a DV-like term in eq. (3.2) is related to the phenomenon of Stokes discontinuities 
[39, 40]. In the case of the ip- and -(//-functions the imaginary axis is a Stokes ray, beyond which, that is in 
the imaginary left-half plane, the formally exponentially suppressed terms combine to yield a potentially 
non-negligible contribution. 
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Figure 1. Integration contour E z ((p) in the complex _B-plane for the FESR (2.15), employing 
zq = 1.5 and the further parameters provided in the text. The part of the contour where the DV 
term (3.3) should be included has been depicted as the solid line, and the remaining one, where the 
asymptotic expansion (3.1) is sufficient, as the dashed line. Some bound state poles as well as the 
continuum cut are also shown as thick dots and line, respectively. 

the -E-plane is displayed in figure 1 for zq = 1.5, together with some bound state poles 
as well as the continuum cut (thick dots and line). The section on which tp'(z) admits an 
asymptotic expansion without DV term is depicted as the dashed line, while the section 
where the DV term (3.3) should be included, is displayed as the solid line. 

Next, figure 2 shows the moment R w (Eq) of eq. (2.15) for the weight function w(E) = 1, 
as a function of zq. As expected, at the locations where E z (ir) crosses a bound-state pole, 
R w (Eq) is discontinuous. Employing the full Adler-type function D{E) of eq. (2.10), it is 
a simple matter to verify that the equality between the phenomenological left-hand side of 
(2.15) in terms of p{E) and the theoretical right-hand side is satisfied. On the other hand, 
when simulating the OPE with the asymptotic expansions (3.1) and (3.2) it is found that 
the DV term (3.3) completely dominates the moment. At the lowest considered zq = 1.2, 
the terms polynomial in 1/z only contribute about 0.3 % and at the highest plotted zq = 5.5, 
this contribution is within the numerical uncertainties of the asymptotic expansion. Hence, 
in the respective case, the moment R w (Eq) is completely saturated by the DV contribution. 

There are two reasons that contribute to this behaviour. On the one hand, the Coulomb 
bound states are zero- width resonances and thus the duality violations are expected to be 
very strong. In such a case the DVs can be considered "maximal". On the other hand, the 
weight function w(E) = 1 does not provide any suppression of the resonance region which 
could soften the strong duality violations. The second issue can be investigated further 
by employing weight functions which provide some suppression of the resonance region, 



-7- 



1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 

Z 



Figure 2. The solid line represents the moment R w (Eq) of eq. (2.15) for the weight function 
w{E) = 1, as a function of Zq. The discontinuities appear at the locations where the contour crosses 
a bound-state pole. The DV contribution (3.3) completely saturates this moment and is visually 
indistinguishable from the solid line. 

for example power- like moments W^{E) = {E — E n ) k {E - Eq). Studying these weights 
for small k = 1, 2, 3, it is found that indeed at low zq (close to the breakdown of the 
asymptotic expansion) the contribution of the DVs is suppressed. However, the dominance 
of DVs again quickly sets in for larger zq, and even at low zq the contribution of DVs is 
still sizeable. The origin of the latter observation can be traced back to the fact that the 
DV term (3.3) penetrates some distance into the complex plane, before the exponential 
decay becomes effective. Thus, even with the weight functions W^ k \E), which nullify the 
contribution on the real axis, the residual contribution from the full contour integration 
can remain sizeable. 

The dependence of the DVs on the resonance structure could be studied further by 
providing the Coulomb bound states with a finite width and varying this width. However, in 
view of practical applications of DV models in analyses of hadronic r decays, the following 
discussion will be continued on the basis of a second model for DVs, already investigated 
in refs. [9, 12, 13]. 

4 A model for light-quark correlators 

Following refs. [9, 12, 13], the structure of DVs shall be investigated in a second model, 
which incorporates constraints from Regge theory [23, 24] on the light meson spectrum, and 
which can be chosen to roughly resemble the physical spectrum of the light-quark vector 
current correlator. The model employed below is a simplified version of the one studied in 
refs. [12, 13], serving all purposes of the analysis discussed in the following. Specifically, it 
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is chosen to take the form 



where 



ny(*) = -4 M ^ (s) )+ const., (4.1) 



and c = 1 "^- (42) 



As compared to refs. [12, 13], the lowest lying vector, that is the would-be p meson, is 
included in the ^-function. The reason for this will be explained below. Furthermore, in 
this work the global normalisation is immaterial and has been dropped. For all remaining 
parameters, in the numerical analysis numbers similar to the ones given in ref. [13] will be 
employed: 

My = 770 MeV, A = 1.2 GeV , a = 0.4 . (4.3) 

In order to arrive at an OPE for the model, we require the expansion of the digamma 
function at large or correspondingly, large \u\. The expansion in question, which is 
however only asymptotic, takes the form [38] 

1 °° B 

ip(z) ~ In z - ^2 2n 2n oo in | arg z\ < it) . (4.4) 

n=l 

In fact, the asymptotic expansion of ip'{z) of eq. (3.1) is of course just equal to the derivative 
of eq. (4.4). Though in principle one would have to further expand eq. (4.4) in terms of 
powers of 1/s, in order to arrive at an OPE-like expansion, this is not necessary for what 
shall be discussed in the following, and thus we stick to the asymptotic expansion in powers 
of 1/z. 

For finite-width resonances, the poles are on an unphysical sheet, and |argz| never 
reaches ir. Still, like in the case of if} f (z), in the region Rez < the asymptotic expansion 
can be substantially improved by making use of the reflection relation [13, 38] 

ip(z) = ip(—z) 7r cot(-7rz) (4.5) 

1 °° B 

z ~°° \nz y ^ 7r [ cot(7rz) ± i] ( Re z < 0, Im z ^ ) . 

2z ^-^ 2n z zn 

71=1 

Again, the logarithm and rational parts of both asymptotic expansions (4.4) and (4.5) agree, 
and it is a simple matter to convince oneself that away from the real axis the additional 
term tt [cot(7rz) dbi] is exponentially suppressed, but takes once more care of the nearby 
poles. 

This is also demonstrated quantitatively in figure 3, where the real (left pane) and 
imaginary (right pane) parts of ip(z) are displayed for z = 1.5 ■ exp(i<£>) as functions of 
93 = arg(z). The solid lines correspond to the full function, while the dashed lines show 
the asymptotic expansion of eq. (4.4), including terms up to n = 5. For positive real z at 
this order the asymptotic series assumes its minimal term. As is evident, the asymptotic 
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arg(z) arg(z) 



Figure 3. Re^z)] (left pane) and Imhy^z)] (right pane), for z = 1.5 • exp(iip), as functions of 
if = arg(z). The solid lines correspond to the full result, while the dashed lines show the asymptotic 
expansion of eq. (4.4), including terms up to n = 5. 

expansion breaks down for <p ^ 2.5. The difference between the total result and the 
asymptotic expansion is to a very good approximation covered by the additional term in 
eq. (4.5), such that taking this contribution into account, the quality of the expansion is 
analogous to the one for ip < tt/2. 

To continue, two finite-energy sum rules for the model (4.1) shall be investigated. For 
an arbitrary, analytic weight function w(s), they take the generic form 

SO 

R w (s ) = [ ^ W (s) pv ( s ) = -L<[^. w ( s )Il v (s), (4.6) 
J s 2tiJ s 

s 

where pv{s) = ImIIy(s + iO)/7r, and the contour on the right-hand side is chosen such that 
it starts and ends at a real so > 0. 

The asymptotic expansion is again studied most cleanly for contours with constant 

absolute value of the argument of the ^-function, \z\ = z$. The corresponding energy 
variable s z ((p) is then given by 

/ M 2 

*z(<p) = -A 2 (W^-^J , (4.7) 

and is displayed in figure 4 for the choice zq = 1.5. Again, the portion of the contour on 
which the DV term should be included is depicted as the solid line, while the remainder, 
on which the expansion (4.4) is admissible, is plotted as the dashed line. In order that the 
contour closes below the cut of IL/(s), that is s 2 (0) < 0, one requires zq > My/ A 2 w 0.41. 
On the other hand this entails that so ^ 1.2GeV 2 . 5 Because we have a model with finite- 

J This is the reason why the p-meson has been included into the ^-function. In the original model of 
ref. [13], the corresponding requirement would have been zo > My /A 2 w 1.6, which would have necessitated 
values for so at least as large as so ^ 4.5 GeV 2 . 
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Figure 4. Integration contour s z (ip) of eq. (4.7) in the complex s-plane for the FESR (4.6), 
employing zq = 1.5 and the further parameters provided in the text. The part of the contour where 
the DV term should be included has been depicted as the solid line, and the remaining one as the 
dashed line. The cut of Ily(s) on the positive real s-axis is also indicated as the thick line. 



width resonances and the poles are located on an unphysical sheet, the angle ipo with 
s z(¥>o) = so is found smaller than ir. For example in the case of figure 4, in which zq = 1.5, 
ifo = 2.972. In the limit of zero-width resonances, corresponding to the large- N c limit, and 
which in the model (4.1) is realised as £ — > 1, (fo would go to tt and the poles would lie 
on the positive real s-axis. On the contrary, for broader resonances which lie further away 
from the real s-axis, the angle <po is smaller, and the contribution of DVs gets reduced. 

The first example of a moment R w (so) corresponding to eq. (4.6) is displayed in figure 5 
for the trivial weight w(s) = 1 and as a function of sq. The oscillatory behaviour of the 
moment results from the presence of resonances, being damped for higher energies. As the 
long- and short-dashed lines, the asymptotic expansion (4.4) is plotted up to the fourth 
and sixth order in z respectively. As is apparent, the asymptotic expansion breaks down 
at about sq ~ 1.5 GeV . Furthermore, large deviations of the asymptotic expansion to the 
exact moment are observed, which are, however, perfectly described by the additional term 
in eq. (4.5). The largest deviation is found at the first minimum around so ~ 1.76 GeV 2 , 
where the DVs amount to about —11% of the total contribution. 

To investigate the influence of weight functions which provide some pinch suppression 
at s = so, in figure 6 the moment R w (so) is displayed for the weight w(s) = (1 — s/so) 2 , 
which has a double zero at s = so, just like the kinematical weight function in hadronic r 
decays. Obviously, now the influence of DVs is much smaller though still clearly visible. 
At so ~ 1.76 GeV 2 , the contribution of DVs amounts to only 1.6% of the total moment, 
and a much stronger suppression is observed towards higher energies. 
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Figure 5. The solid line shows the moment R w (sq) of eq. (4.6) for the weight w(s) = 1 and plotted 
as a function of sq. The long- and short-dashed lines correspond to the power- like asymptotic 
expansion up to fourth and sixth order respectively. 
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Figure 6. The solid line shows the moment R w (so) of eq. (4.6) for the weight w(s) = (1 — s/sq) 2 
and plotted as a function of sq. The long- and short-dashed lines correspond to the power- like 
asymptotic expansion up to fourth and sixth order respectively. 

5 Conclusions 

In phenomenological analyses of QCD, the operator product expansion is a widely used tool. 
While at sufficiently high energies the legitimacy of its application is unquestionable in the 
Euclidian domain, the continuation of the OPE towards the physical, Minkowskian region 
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is inflicted with the appearance of duality violations. The presence of DV contributions, 
besides the OPE, is due to the physical bound states or resonances, whose physics cannot 
be captured by the OPE alone. 

In the present article, the structure and impact of DVs has been studied on the basis of 
two models. Firstly, the quantum mechanical Coulomb system was investigated, which can 
serve as a model for heavy quarkonia, and which has not been considered in the context of 
finite-energy sum rules and DVs before. Secondly, a simplified version of a model already 
discussed extensively in refs. [9, 12, 13] was used for comparison with the features found in 
the Coulomb model. The latter model incorporates constraints from large- iV c QCD, as well 
as Regge theory, and captures the central features of the light-quark vector or axialvector 
correlators. 

Qualitatively, the physical picture behind the appearance of duality violations can 
be summarised as follows: in the Euclidean region, all effects of resonances or bound- 
state poles have been sufficiently smoothed out, quark-hadron duality is expected to work 
well, and the OPE provides a satisfactory description of current correlation functions. 
When analytically continuing through the complex energy plane towards the physical, 
Minkowskian region, the presence of resonance poles becomes more and more prominent, 
reflecting itself in a breakdown of the OPE. Mathematically, the DV contributions show 
up as formally exponentially suppressed terms in an asymptotic expansion, which are, 
however, enhanced by the nearby poles, and thus can become numerically relevant. 

The strength of effects resulting from DVs can be considered "maximal" if the bound 
states have zero width and the poles lie on the Minkowskian axis. This case is for example 
realised in the Coulomb model or when N c in eq. (4.2) of the second model is taken to 
infinity [12]. On the other hand, the influence of DVs is suppressed for resonances with 
a finite width and becomes smaller when the resonances are broader and the poles lie 
further away from the physical, Minkowskian axis. This is also seen from the oscillations 
in the physical spectrum whose features cannot be described adequately by the OPE. These 
oscillations are more prominent for narrower resonances while they become weaker if the 
resonances get broader. 

If the physical resonances are sufficiently broad, and thus lie sufficiently far away from 
the physical axis, a simplified description of DVs as an exponentially damped oscillatory 
function appears admissible. Assuming s to be large enough for the OPE to make sense, 
motivated by the ^-function model such a simplified Ansatz was advocated in refs. [13, 14] 
for the duality violating spectral function Pdv( s ) i n the form 



where k, 7, a and /3 are a priori free parameters which in a specific channel have to be 
extracted from experiment, and which bear no immediate QCD interpretation. The DV 
contribution to a moment is then assumed to be given by 




(5.1) 
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In practical applications of QCD finite-energy sum rules, however, the DVs naturally 
are required on the complex contour and thus a closely related Ansatz, which corresponds 
to the continuation of eq. (5.1) to the right complex half-plane, may also be considered: 

n DV (s) = TTKe-~< s±i{a+ P s) (Res > 0, Ims> 0), (5.3) 

with 7 and (3 positive. Of course, the IhiIIdv( s ) ° n the positive, real s-axis reproduces the 
spectrum of eq. (5.1), but the Ansatz also reflects the form of the exponential suppression of 
the DV term given in eq. (3.3). Both Ansatze (5.1) and (5.3) would be totally equivalent, 
if the DV contribution would be active on all of the complex contour. However, since 
the DV term should only be included on part of the contour (the solid-line sections in 
figures 1 and 4, or, more generically, on the right-half s-plane Res > 0), both choices 
differ by an exponentially suppressed contribution. As a matter of principle this should 
be numerically irrelevant, but it could play a certain role if in practical applications the 
exponential suppression is not very pronounced. This will be studied in the future in 
investigations like the ones of refs. [14, 41]. 

To conclude, in the presented article observations about duality violations appearing 
in asymptotic expansions were summarised on the basis of two models. One of them has 
already been discussed in the literature before, while the other, the Coulomb system, has 
not been studied with respect to DVs previously. The qualitative findings of this work 
should prove useful for phenomenological QCD studies in the future. 
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